UNSTEADY OUTFLOW MODELS FOR COSMOLOGICAL GAMMA-RAY BURSTS 
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Abstract 



The 'event' that triggers a gamma ray burst cannot last for more than a few seconds. 
This is, however, long compared with the dynamical timescale of a compact stellar-mass 
object (~ 10 -3 seconds). Energy is assumed to be released as an outflow with high mean 
lorentz factor V. But a compact stellar-mass collapse or merger is, realistically, likely to 
generate a mass (or energy) flux that is unsteady on some timescales in the range 10 -3 
- 10 seconds. If T fluctuates by a factor of ~ 2 around its mean value, relative motions 
within the outflowing material will themselves (in the comoving frame) be relativistic, and 
can give rise to internal shocks. For T ~ 10 2 , the resultant dissipation occurs outside 



the 'photosphere' and can convert a substantial fraction of the overall outflow energy into 
non-thermal radiation. This suggests a mechanism for cosmological bursts that demands 



2 ■ less extreme assumptions (in respect of T-values, freedom from baryonic contamination 
. etc) than earlier proposals. 



1. Introduction 

Gamma ray bursts are clearly the 'signal' of an energetic event lasting (typically) no more 
than a few seconds. The energy output in the gamma rays themselves, if these events are 
at 'cosmological' distances, would be up to 10 51 ergs. Each burst would then, probably, 
involve the collapse of a stellar-mass object, or the coalescence of a compact binary (see 
Hartmann, 1993, or Paczyhski , 1994, for an overall review). The main theoretical challenge 
is to understand how the energy can be converted, with adequate efficiency, into gamma 
rays with a non-thermal spectrum 

A widely recognized problem is that if the rest mass energy of entrained baryons 
exceeded even 10" 5 of the total energy, the associated opacity would trap the radiation so 
that it was degraded by adiabatic expansion (and thermalized) before escape, e.g. Cavallo 
and Rees, 1978, Paczyhski , 1990. The energy would therefore be transformed into kinetic 
energy of bulk relativistic outflow. This problem arises if the 'event' is approximated as 
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an instantaneous fireball, or as an outflowing wind which is 'steady' over the entire burst 
duration. In earlier papers, we have discussed how kinetic energy can be reconverted into 
gamma rays by relativistic shocks which form when the ejecta run into external matter 
(either ambient interstellar matter, or a non-relativistic outflow preceding the 'event' it- 
self). We showed that acceptable models required bulk Lorentz factors V of order 1000 - 
still high, but allowing much more baryonic loading and opacity than would be tolerable if 
the energy escaped directly from a 'simple' fireball before adiabatic losses had attenuated 
and thermalized it. 

In this paper we show that the constraints are eased still further if we adopt a 
less idealized picture of the 'event' itself. The production of the relativistic ejecta will be 
spread over a finite duration; moreover the physical conditions determining rj = E/Mc 2 
(or 7] = L/Mc 2 , the ratio of radiation and magnetic energy to rest mass, which gives an 
upper limit for V) will not be steady throughout the event. Fast (higher if) ejecta can catch 
up with slower material that was ejected earlier: kinetic energy can then be reconverted 
into energetic particles (and thence into gamma-rays). This suggests a mechanism for 
'cosmologicaP bursts that can operate for much lower values of rj (i.e. higher loading 
factors) than were previously believed necessary. 

In these more realistic models, dissipation happens whenever internal shocks de- 
velop in the ejecta - it need not await the deceleration by sweeping up ambient external 
matter. Simple relativistic kinematics show, however, that the internal shocks do not de- 
velop until the ejecta have attained sufficiently large distances that the resultant radiation 
can escape without thermalization or adiabatic losses. 

We show how variations in the terminal speed, resulting from inhomogeneous or 
time- varying conditions around the central object, can yield efficient production of gamma 
rays, and generally a complex time-structure. We outline the general kinematics in §2. We 
then discuss in §3 the nature of the dissipation, and the distinctive role of strong fields 
in 'magnetically-dominated' outflows. The relation of these ideas to gamma-ray burst 
phenomenology is outlined in §4. 

2. Kinematics in an Unsteady Relativistic Outflow 

2.1 An Illustrative Example 

We postulate an outflow persisting (typically) for a few seconds. But instead of assuming 
this to be a steady wind we suppose that it is irregular on much shorter timescales. This 
may be a 'wind' from a high-S newly formed pulsar (Usov 1992,94), or the debris (again 
probably highly magnetized) flung off from a compact binary during the complex dynamics 
of its coalescence into a black hole and surrounding disc (Narayan, et al. , 1992, Davies, 
et al. , 1994). 

Whenever part of the ejecta 'catches up' with other material ejected earlier at a 
lower Lorentz factor, an internal shock forms, which dissipates the relative kinetic energy. 
To illustrate the basic idea, suppose that two 'blobs' of equal rest mass, but with different 
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Lorentz factors Yi and T 2 (with r 2 > > 1) are ejected at times t\ and t 2 , where 
t<i — t\ = t var . In the case of highly relativistic ejecta, the shock develops after a distance 
of order ct var YiY2- For high Lorentz factors, therefore, the shock takes a long time to 
develop, even if V 1. This is, of course, because the distance that must be caught up 
is (in the 'laboratory' frame) of order ct var , but the speeds all differ from c by less than 

i/r?. 

For example, suppose that the Lorentz factor of the outflow is, on average, 100, 
but varies from 50 to 200 on a timescale t var . The velocity differences are of order 10 _4 c, 
so the distance for the shock to develop in the lab frame is 10 4 ct„ ar . The reconversion of 
bulk energy can nevertheless be very efficient: when the two blobs share their momentum, 
they move with Y fi na i = v^FiF^ , so the fraction of the energy dissipated is 

e = (y 1 + r 2 - 2 v / f^)/(r 1 + r 2 ) . (i) 

For the previous numerical example, the efficiency would be 20%. High efficiency does 
not, therefore, require an impact on matter at rest; all that is needed is that the relative 
motions in the comoving frame be relativistic -i.e. Y 2 /Yi > 2 (c.f. Rees, 1978, for an 
application of this argument in a different context). 



2.2 An Unsteady Wind 

Suppose that the mean outflow (over the few seconds that a typical event lasts) can be 
characterized by a steady wind with given mean values of L w and rj = L w /Mc 2 . We then 
envisage that the value of rj (or L w ) is unsteady. 

The mean properties of the wind determine the average bulk Lorentz factor 

r ~ / ( r / rz ) ' for r ^ Ts > (2) 

\ i] , for r ^ r s . 

Here ri is the size at the base of the wind (in 'young pulsar' models, c.f. Usov, 1992, this 
would equal cP/2ir, where P is the period). The Lorentz factor saturates to Y ~ rj at a 
saturation radius r s jr\ ~ rj where the wind energy density, in radiation or in magnetic 
fields, drops below the baryon rest mass density in the comoving frame. For a given 
baryonic mass loss M the photospheric radius where the wind becomes optically thin to 
Thomson scattering is 

r ph = Mk/(4tvcY 2 ) = 1.2 x 10 12 L 51 r/^ 3 cm , (3) 

where r/2 = (r//10 2 ) and L 5 i = (L^/10 51 ergs). The above equation holds provided that rj 
is low enough that the wind has already reached its 'terminal' Lorentz factor at r p h- This 

requires rj ^ rj m ~ 10 2 L§{ 4 t V ar ' 4 , if one takes r\ ~ ct var . 

When a compact binary coalesces, the ejecta would obviously be very messy. The 
characteristic timescale at the 'base' of the wind is of order td yn = 10~ 3 s. (The breakup 
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angular speed of a pulsar is of the same order) . Large-amplitude variations could occur on 
any timescale longer than this. (Even more rapid variations are also possible, especially if 
there is a strong tangled magnetic field in the wind, e.g. §3). 

If the value of 77 at the base increases by a factor ^ 2 over a timescale t var , then 
the later ejecta will catch up and dissipate a significant fraction of their energy at some 
radius > r s given by 

rd ~ ct var r] 2 ~ 3 x 10 14 t var rj2 cm , (4) 

where t var is in seconds. Dissipation, to be most effective, must occur when the wind is op- 
tically thin - otherwise it will suffer adiabatic cooling before escaping, and be thermalized. 
(Outside r s and r p h, where radiation has decoupled from the plasma, the sound speed will 
be far below c; this also guarantees that the relativistic internal motions in the comoving 
frame will lead to shocks in the gas). This implies the following lower limit on 77: 

77 ^ 3 x lO^fC/ 5 . (5) 

This simple estimate doesn't take into account the extra pairs that may result from dissi- 
pation. We comment further on this in §3. 

The physical conditions in these shocks qualitatively resemble those in the reverse 
shocks behind the blast waves discussed by Meszaros and Rees (1993); however the densities 
and magnetic fields are higher, and the coling more efficient, because r is smaller. (The 
internal Lorentz factors are modest, because, to first approximation, all the material is 
comoving outwards at more or less the same speed) . 

There will be variations with a range of t var . Rapid fluctuations are dissipated at 
smaller r than those on longer timescales. There would therefore be a dependence of the 
spectrum on the characteristic variability timescale . The radiation processes depend on 
the magnetic field strength. We discuss this next, because the field may also be dynamically 
important. 

3. The Role of the Magnetic Field. 

Ultra-intense magnetic fields are expected either in 'young pulsar' models (Usov 1992) or if 
the field builds up towards equipartition by differential motions (Narayan et al 1992) or a 
convective dynamo (Duncan and Thompson, 1992). Collapsing or coalescing neutron stars 
may generate fields as high as Bi ~ 10 16 G. Magnetic stresses could, indeed, be dynamically 
dominant over the radiation: the ratio of magnetic energy to rest-mass energy at the base 
of the wind (r = 77) would then determine the effective value of 77. Even a magnetic field 
that was not dynamically-dominant would still be important in ensuring effective cooling 
by cyclotron and synchrotron emission. 
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If the Poynting flux provides a fraction a of the total luminosity L at the base 
of the wind (at r/ ~ ct var ) the magnetic field there is B[ ~ 10 10 a 1 1 ' 2 L\'^t~^ r G. The 
comoving magnetic field at the dissipation radius (4) (which always is outside r s ) is 

B d = BM/rsfirJra) ~ lOV/ 2 ^ 2 *^" 3 G . (6) 

If the electrons are accelerated in the dissipation shocks to a Lorentz factor 7 = 10 3 73 
the ratio of the synchrotron cooling time to the dynamic expansion time in the comoving 
frame is 

{t ay /tex)d ~ 5 x 1(T 3 a' 1 L^^H var ril , (7) 
so a very high radiative efficiency is ensured even for t var as high as seconds. 

If the shock dissipation leads to photons whose energy in the comoving frame 
exceeds 1 Mev, then there is the possibility of extra pair production from photon collisions; 
this would modify eq.(5). The pairs could increase the effective photospheric radius by a 
factor x = (m p /m e ) x (radiative efficiency) x (fraction of radiation going into photons above 
1 Mev). For the dissipation to occur outside any possible pair-dominated photosphere (a 
requirement that may actually be unnecessary if the pairs annihilate on a shock cooling 
timescale) 77 would need to be higher than in eq.(5) by a factor x 1 ^. So even if x were 
(say) 100, the minimum required rj would go up by by less than a factor of 3. 

It is clear from eq.(7) that a magnetic field can ensure efficient cooling even if it is 
not strong enough to be dynamically significant (i.e. even for a C 1). If, however, the field 
is dynamically significant in the wind, then its stresses will certainly dominate the (pre- 
shock) gas pressure. Indeed, in a wind with a = 1 the magnetosonic and Alfven speeds may 
remain marginally relativistic even beyond r s if the field becomes predominantly transverse. 
In this extreme case, magnetic fields could inhibit shock formation unless 77 varied by much 
more than a factor of 2. On the other hand, the presence of a dynamically-significant and 
non-uniform field could actually drive internal motions leading to dissipation even in a 
constant-77 wind. Except in the special case of an aligned dipole, the magnetic field would 
have reversals on a scale of order r\ ~ cP/2n = 5 x 10 6 P_3 cm. If the field inside r\ has a 
complex (non-dipole) structure, the reversal could be even smaller. (Thompson, 1994, has 
discussed a detailed model where the resultant Alfven waves are dissipated via Compton 
drag when the scattering optical depth is still large.) 

4. Phenomenology and Discussion 

An unsteady (and probably magnetized) wind or fireball has the advantage that it can 
accomodate a larger amount of baryon contamination than in previous models, while still 
producing a nonthermal gamma-ray burst via synchro-Compton radiation (as in the reverse 
shock of Meszaros , Laguna and Rees, 1993; c.f. also Katz, 1994) from electrons accelerated 
in the shock dissipation region beyond the photosphere. 

If the energy were released as an 'instantaneous' fireball, or as a steady wind, 
and transformed into kinetic energy by adiabatic expansion, then efficient reconversion 
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into gamma rays occurs when (but not until ) enough external matter has been swept up 
to decelerate it. In earlier papers, we have shown how the resultant shocks could give 
rise to gamma-ray bursts provided that the Lorentz factors are of order 1000; complex 
time-structure must then be put down to irregularities in the ambient medium. 

In the present paper, we have assumed (undoubtedly more realistically) that the 
energy release is complex and irregular, and shown that this assumption admits the extra 
possibility that 'internal shocks' can dissipate a substantial fraction of the kinetic energy 
before the ejecta encounter the ambient medium. The compact object triggering the burst 
is likely to have a characteristic dynamical timescale of only about td yn ~ 10 -3 s (of order 
tvar ~ ti/c in §2). But the energy release is likely to be more prolonged, determined 
by, e.g., magnetic spindown, disk viscosity or neutrino diffusion timescales, depending on 
the model. The energy flux in the outflow, or the value of 77, could then fluctuate on all 
timescales from tdyn up to the overall duration of the energy release, which could be many 
seconds. (Indeed, some models - e.g. rapidly spinning pulsars with complex non-dipole 
fields - permit irregularities on timescales even shorter than td yn )- 

In comparision with shocks involving external ambient matter, the internal shocks 
form at smaller radii, and in regions of higher density. The Lorentz factors (and values of 
77) needed in order to get efficient dissipation, and short 'observer frame' timescales, can 
then be somewhat lower. However, values of r/ <S 30(L 51 /t t;ar ) 1 / 5 (eq.[5j) are still required 
in order to ensure that the shocks form outside the photosphere. 

This implies wind mass losses M ^ 3 x 10" 2 Lc~ 2 ~ 10" 5 L 51 M s _1 . The total 
(isotropic) mass loss expected from an unmagnetized collapsing core or merging compact 
binary is of order 1O -3 M0 s _1 (e.g. Meszaros and Rees, 1992, Woosley, 1993). However, 
along the rotational (or binary) axis centrifugal forces may significantly reduce the baryon 
losses, while in a strongly magnetized object mass loss is expected only from the open polar 
field lines, representing a fraction 10 -2 of the total area, so M ^ 1O _5 M0 is reasonable. 
(The value of M implied by our required value of i] (~ 10 2 ) are high enough to ensure 
that the MHD approximation applies to any magnetic field in the wind; this contrasts with 
Usov's (1992, 1994) proposal, which would require a very much lower M). 

The observed bursts are remarkable for their disparate and complex time- 
structure, and this may be a consequence of how the development of internal shocks is 
subject to irregularities in the outflow. 

In summary, the implications of this work are 

(i) The short timescales (and adequate efficiencies) do not need such high values of 77 as 
our earlier blast wave models (which were themselves much less demanding in this respect 
than pre-1992 'fireball' models). 

(ii) The time structure could be complex, being dependent on the time history of the 
lorentz factor. The dissipation associated with shorter timescales would tend to occur at 
smaller radii. 

(iii) Since the observed gamma rays from each part of the wind are concentrated, owing to 
aberration, into an angle of order 1/r/, our discussion can be straightforwardly extended 
to a 'beamed' or jet-like geometry. The broad range of burst morphologies could then , at 
least in part, be due to the axis and the boundaries of jets having different Lorentz factors 



6 



or internal variability. 

(iv) The energy sources in a cosmological context could be either stellar collapse (Usov, 
1992, Woosley, 1993) or compact binary mergers (Paczyhski , 1986, Eichler, et al. , 1989, 
Narayan, Paczyhski and Piran, 1992). 
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Figure Caption 

Fig. 1: Wind regimes as a function of n = L/Mc 2 and r/ctd yn - The lines above 
which the wind becomes optically thin to scattering (due to M) has a slope -1/3 or -1/2 
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for L or M constant after saturation, and is fixed before saturation. The lines of t var 
proportional to td yn have a slope 1/2. GRBs from self-consistent unsteady winds are in 
the triangular region lying below the line t var ~ ^dyni and above the line t — 1 representing 
the photosphere. 
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